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ABSTRACT 

We used the Magellan adaptive optics (MagAO) system and its VisAO CCD 
camera to image the young low mass brown dwarf companion CT Chamaeleontis 
B for the hrst time at visible wavelengths. We detect it at r', i', z', and Ys- With 
our new photometry and Tg// ~ 2500 K derived from the shape its A—band 
spectrum, we hnd that CT Cha B has Ay = 3.4 ± 1.1 mag, and a mass of 14 — 24 
Mj according to the DUSTY evolutionary tracks and its 1 — 5 Myr age. The 
overluminosity of our r' detection indicates that the companion has signihcant Ha 
emission and a mass accretion rate ~ 6 x lO“^°M 0 /yr, similar to some substellar 
companions. Proper motion analysis shows that another point source within 2" of 
CT Cha A is not physical. This paper demonstrates how visible wavelength AO 
photometry (r', i', z', Ys) allows for a better estimate of extinction, luminosity, 
and mass accretion rate of young substellar companions. 

Subject headings: brown dwarfs — planetary systems — stars: individual (CT 
Cha) — stars: pre-main sequence — instrumentation: adaptive optics — planets 
and satellites: individual (CT Cha B) 
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INTRODUCTION 


As more and more brown dwarfs and planetary companions are being discovered, char¬ 
acterizing them in the visible regime yields a more complete pictnre of the spectral energy 
distribntion (SED) and more insight into physical properties as well as formation scenar¬ 
ios. For instance, a better estimate of extinction helps to derive bolometric Inminosity and 
mass—especially for yonng objects (< 10 Myr) which may still have primeval dnst and 
gas aronnd them and snffer significant obscnration. However, extinction is problematic to 
measnre becanse most of the high-contrast adaptive optics (AO) observations are done in 
the near-infrared, which is ~ 10 times less sensitive to dust at K versus V. One simple 
treat ment is to assume th at the companion has the same amount of extinction as its host 
star fjPatience et al.ll2012[) . since in the early stages of star formation the binary might be 
embedded in a common envelope. For more evolved, fragmented systems both components 
may have their own disks, so this assumption might be invalid. Ideally one would like 
to acquire visible spectra or at least broad-band visible photometry to supplement near- 
IR measurements because visible wavelengths are a better probe for dust extinction. Yet 
high contrast opti cal observations on companions are very rare due to decreased contrast 
(iMales et al.l 1201 411 and the difficulty of correcting atmospheric turbulence at visible (defined 
here as A < l/im) wavelengths. We therefore need an advanced AO system which can work 
in the visible to suppress the halo. 

Here we present the hrst optical AO photometry of CT Chamaeleontis system with 
the Magellan adaptive optics (MagAO) syste m, a powerful new 585 - element AO system 


commissioned on the 6.5 m Cla y Telescope flClose et al.l I2013L I2014J: iFollette et al 


Wu et ah 2013 : Males et al.ll2014 h CT Cha A, a K7 classical T Tauri star f Weintraub 


2013 


1990 


Gauvin fc Stromlll992h. is located in t he Chamaeleon I star-forming regio n. This region is 


close (~ 160 pc; IWhittet et al.l 119971: iBertout et al.l Il999l: iLuhmanl 1200811 and is as young 
(median age ~ 2 Myr; If jnhma.nl 1200411 as the Taurus sta r-forming region and I C 348. It 


also has relatively low extinction (typical Ay < 5 mag; ICambresv et al.lll997l) . enabling 
clear view of youn g stars. The companio n CT Cha B at 2.67"(430 AU) projected separation 
was discovered by Schmidt et ah f 2008[l in their VLT NACO survey. Based on its near- 
IR spectrum, the companion was estimated to be an M8-to-L 0 (Tpff ~ 2600 K) low mass 
(~ 17 Mj) brown dwarf with Ay ~ 5.2 mag. ISchmidt et al.l fl2008l) also imaged another 


*This paper includes data gathered with the 6.5 m Magellan Clay Telescope at Las Campanas Observatory, 
Chile. 
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closer object termed “cc2”, whose true nature has remained puzzling fISchmidt et al.l 12009 


Robberto et al.ll2012[) . In this paper we present new visible AO observations providing a 


better measurement of Ay and of the mass of CT Cha B. Our accurate astrometry allows 
us to determine that cc 2 is, in fact, a background source. 


OBSERVATIONS and REDUCTION 


MagAO observations with its VisAO camera fj Close et ahl l2013l: Males et al.l l2014l) at 
r' (0.62/im), i' (0.77/im), z' (0.91/rm), and Vs (0.98/rm) were performed on 2013 April 6 (UT) 
during the second commissioning run. Seeing was stable, ranging from 0.6"to 0.8". We 
locked the AO system on CT Cha A (R ~ 12 mag) at 100 modes and 625 Hz@ The achieved 
FWHMs were 0.1", 0.08", 0.06", 0.06" for r', i', z', Ys, respectively. Strehl ratios were 
low due to only correcting 100 modes, since the guide star was somewhat faint for VisAO. 
We obtained saturated images to boost signal-to-noise ratio, with unsaturated datasets for 
relative photometry (top row in Figured]). As a young accreting star, CT Cha A varies its 
brightness by ^ 0.3 mag in the nea r-IR and ~ 1 mag in the optical flBatalha et al.l Il998 


Chez et al.lll997l: iLawson et al.lll996h . In order to calibrate the brightness of CT Cha A, we 


also obtained absolute photometry by observing the optical photometric standard star LTT 
3864. 

We carried out standard data reduction with IRAF tasks. After dark-subtraction, flat- 
held correction, and cross-correlation between frames, we rotated counterclockwise the satu¬ 
rated data by 89.11°-|-parallactic angle to make north up and east left (middle row in Figure 
d]). This was following the cal ibration of VisAO camera based on astrometry of the Trapez¬ 
ium cluster fjClose et al.l 120131) . We further rotated each of these images by 20, 40, ..., 340° 
and took the median of them to approximate the halo of the primary. Then we subtracted 
the halo from the original images to further bring out any faint point source object (bottom 
row in Figured]) without any loss of hux from self-subtraction. Anisoplanatic ehects are still 
small at this small separation (< 3"). Only our bluest hlter (r') showed some sign of aniso- 
planatism, so we smoothed the reduced r' image with a gaussian (width = 2 x PSF FWHM) 
to enhance signal-to-noise ratio. Next we constructed a master PSF from unsaturated im¬ 
ages of A, used it to £t CT Cha B prohle in the deep images, and measured the PSF htting 
flux with the DAOPHOT allstar task. We note here that the on-axis CT Cha A PSF was 
still an excellent £t to the 2.7" off-axis PSF of B. Table dl summarizes our observations and 


§The faintness of this guide star prevented us from using all 378 modes at 1000 Hz which typically requires 
i? < 10 mag guide stars. 
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PSF fitting photometry on CT Cha B. Uncertainties of B are <0.1 mag for F, z', and Ys, 
and ~ 0.2 mag for r' due to lo w S/N. Near-IR (J, H , Ks) photometry and the ~ 0.3 mag 
uncertainty were adopted from ISchmidt et al.l (120081) . 


To increase the accuracy of our astrometry, we also corrected for image distortion (< 
15 mas). The exact forr nulae to co r rect a ny residual distortions for separation {6x,6y) 
from {X,Y) are listed in IClose et al.l (1201311 and reproduced here: truefo = measured^ — 
6dx X [measured^ 1/110.0 and true^^^ = measured^^^ — 6dy x |measured5j^|/44.5, where 6dx = 
-0.00038921676 x (X - 512) + 0.00084322443 x (U - 512) and Sdy = -0.00025760395 x 
(X — 512) — 0.0024045175 x {Y — 512). We also retrieved archive HST and VLT NACO data 
for proper motion analysis. The HST data were already re d uced by the OPUS pipeline and 
the MultiDrizzle software, as described in iRobberto et al.l (120121) . We reduced NACO raw 
frames by shift-and-add, without flat-held correction and dar k-subtraction. NACO image 
distortion is very small and only up to 3 mas at held edges (jNeuhauser et ah 2008). The 
error budget of our measurements includes platescale and centroid uncertainties. Table [2] 
lists our astrometric measurements. 


3. RESULTS 


3.1. Optical Images 


Figured] shows the CT Cha system seen with our broad-band hlters, with CT Cha B and 
cc2 visible in all four. This is the hrst optical detect ion of CT Cha B, as it was not detected 
in previous HST narrow-band optical observations (IRobberto et al.ll2012[) . Judging from its 
color, cc2 is rel a tively blue (r' — i' = 0.9) so unlikely to be another low-mass companion. 
Robberto et al.l (120121 1 also speculated that a faint “object” seen at [Of] 1.5" to the south 
of CT Cha A could be real, but we cannot conhrm any other faint object in our images, 
especially at r' where a narrow band [Of] or Ha source might have been visible. Thus, it is 
unlikely to be a real object. 


3.2. Astrometry 


The nature of cc2 is not fully settled in literature. ISchmidt et al.l (l2009l) presented 


two-ye ar astrometry, showing that it is likely to be a background star. But IRobberto et ah 


(I 2 OI 2 I) suggested that cc2 may be physically associated based on their single epoch HST 
observations. We measured the positions of cc2 and CT Cha B in images taken by various 
instruments over ~ 7 year time span (Figure |2|). Signihcant common-proper motion has been 
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found for CT Cha B, confirming it is physically bound. However, we detected a significant 
non-common ~ 15.4 mas/yr northwestward motion for cc2, unambiguously demonstrating 
that it is not a co-moving companion but instead a background star, and not a member of 
Chamaeleon I. 


3.3. SED Fitting and Derived Properties 

3.3.1. Effective Temperature 


To further narrow down the u ncertainty of Tp ff. w e retrieved the spectrum taken with 
the VLT SINFONI spectro graph (Schmidt et ah 2008 1. and calculated the H 2 O-K 2 index 
following the prescription of iRoias-Avala et al.l (120 12[) . Assuming solar metallicity, we found 
that for CT Cha B this index is almost independent of extinction, ranging from 0.65 to 0.66 
for Ay = 0 to 5.5 mag. In Figure [3l we plotted the variation o f the index with a range of 
surface gravity for 2000 — 2800 K BT-Settl atmospheric models (jAllard et ah 2011 1. Within 
this temperature range, H 2 O-K 2 index is rather insensitive to log g. Hence we are free to use 
it for young cool objects like CT Cha B. Our best £t to the index corresponds to a spectral 
type M9 ± 1 with T^ff = 2500 ± 100 K. 


3.3.2. Extinction, Bolometric Luminosity, and Mass 


CT Cha B was previously estimated from near-IR sp ectroscopy to have a n extinction 
higher than its host star {Ay =5.2 mag versus 1.3 mag; ISchmidt et al.N2008fl . Our data 
benefit from the fact that visible wavelengths are more sensitive to dust extinction, so we 
can determine Ay at higher precision with MagAO’s VisAO camera. 


We applied multiple values of Ay following the extinction law in iFitzpatrickI fll999l) to 
redden the 2500 K BT-Settl synthetic spectra normalized at Ks (Figure 0]). Minimization 
of the reduced is based on the reddened models £t to the observed i', z', Ys, J, and H 
photometry (black points in Figure |5]). We found that while the result is independent of 
surface gravity, Xr remains high even after including the ±0.3 mag uncertainty at Ks. Some 
systematic errors may come into play. For example, the adopted extinction law might be 
invalid due to grain growth in the disk, or there could be multiple dust components. On the 
other hand, scattered lig ht from the disk or outflow gas may contribute to our i' photometry, 
as in the case of R Mon fjClose et al.lll997() . In this picture, blue light follows indirect paths 
to the observer, avoiding passing through the disk and making our extinction estimate likely 
a lower limit. Another possible cause for higher Xr is fh^t an overall offset ~ 0.5 mag might 
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exist between the visible and near-IR data because they were taken on different nights and 
CT Cha A is a well-known variable. Finally, the companion itself could also be variable 
in the visible just like the primary due to accretion. In any case, with no prior knowledge 
of the material around CT Cha B, our current data yield a best fit to a lower extinction 
Ay = 3.4±1.1 mag. We plotted the reddened synthetic spectrum together with the observed 
photometry in Figured 


We followed the approach in iHillenbrandl fjl997|) to calculate the bolometric luminos¬ 
ity. We converted our i' photometry to Cousins I r , de-r edde ned it by Ay = 3.4 mag, ap¬ 
plied the bolometric correction from iTinnev et al.l fll993l) and iBesselll fjl995[) . and obtained 
Iok(L hni/ Li^) = —2.68 ± 0.2 5. As a comparison, we also calculated L^oi from the K flux 
following IClose et al.l fj2007|) and had a similar value \o?^{ LhnilLcA = ~2 - 71 ± 0 .20. Both 
values are consistent with \og{Lhoi/L q) = —2.68 ± 0.21 in I Schmidt et ahl (120081) . We also 
calculated CT Cha B’s radius u sing L oc and obtained ~ 2.4Rj. Then we applied 

the DUSTY evolutionary tracks (IChabrier et al.ll2000l: iBaraffe et al.ll200ll ) to derive a mass 
estimate of ~ 14 — 2AMj based on the ~ 1 — 5 Myr age and Tg// (Figure [6]). Therefore, CT 
Cha B is most likely a very low mass brown dwarf, just above the planetary mass limit. 


3.3.3. Accretion Rate 


Pa—d emission, an accretion signature, has been seen in CT Cha B’s J—band spectrum 
fjSchmidt et al.ll2008l ). Since CT Cha B is widely separated from the host star, it may harbor 
its own disk and still be actively accreting at this time. Figure |5] shows that our r' detection 
is about 20 times brighter than its predicted continuum. This significant r' excess seems to 
imply strong Ha emission from accretion, allowing us to calculate the mass accretion rate . 
Attributing > 95% of the r' flux to Ha and following the prescription of IClose et al.l (120141) . 
we estimated M ~ 6 x 10“^°MQ/yr, which is reasonable as it implies that a few Mj of 
brown dwarf mass could be accreted in a few million years at the end of the gas-rich disk 
phase. The accre tion rate we derived is also consistent with recent HST observations by 
Zhou et al.l (120141 ). who measured M ~ 10“^^ — 10“®MQ/yr for three substellar companions 


CSC 6214-210 B, GQ Lup B, and DH Tan b based on their optical excess. 


3.4. Implications 

The different extinction between the primary and the secondary may imply that both 
objects have their own disks likely with different inclination angles, resembling conceptually 
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the configuration of HK Tan A and B flJensen fc Akeson 


20141) . CT Cha B’s r' excess, 


together with other accreting objects in IZhou et al.l fl2014i) . suggest that accretion disks 
could be common among young low-mass companions and favor the “star-like” formation 
via gravitational collapse and fragmentation of molecular clouds. The survival of these 
significant disks also implies th at substellar comp anions form near their current locations 
rather than being ejected there flKraus et al.ll2014j) . Strategic Ha surveys such as MagAO’s 
ongoing Giant Accreting Proto-planets Survey (GAPplanetS) may have the potential to 
probe ^ lAfj accr eting giant planets and shed light on the earliest stage of planet formation 
I Glose et ah 2014 1. 


4. Summary 

MagAO observations on GT Gha at r', i', z'^ and Ys have improved the accuracy of the 
extinction towards GT Gha B. The companion is detected in all of our optical filters, whereas 
no detections were made by HST. It is over-luminous at r', indicating active accretion at a 
rate of M ~ 6 x The H 2 O-K 2 index derived from the Ks spectrum is consistent 

with a Te// = 2400 — 2600 K brown dwarf. Using the BT-Settl model, we show that GT Gha 
B is best fit by Ay = 3.4 ± 1.1 mag, which is lower than previous estimate and translates 
to a mass estimate of 14 — 24 Mj based on the DUSTY tracks. We do not see the faint 
southern [01] source seen in previous HST observations, so it is unlikely to be real. Finally, 
our astrometry on cc2 is incompatible with a previous claim that it is a co-moving object. 
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Table 1. MagAO Photometry on CT Cha B 


Filter 

tsat 

(sx#) 

iunsat 

(sx#) 

m 

(mag) 

Am’’ 

(mag) 

Fx^ 

(10“^^ erg s“^cm“^/im“^) 

r' 

20 X 75 

2.27 X 13 

21.71 ±0.21 

9.80 

0.51 ±0.10 

i' 

30 X 52 

2.27 X 14 

20.32 ±0.09 

9.22 

1.00 ±0.09 

z' 

15 X 13 

2.27 X 13 

18.46 ±0.08 

7.96 

3.46 ±0.26 

Ys 


30 X 35 

17.94 ±0.09 

7.59 

4.55 ±0.41 


'^Relative to CT Cha A, measured with PSF htting photometry. 

’’Calibrated with the standard LTT 3864 taken at a similar air mass to CT Cha in 
photometric conditions (errors are < 1% due to different air masse s fro m model at¬ 


mosph ere). We used photometry calibrations of IClose et al.l (120131) and Males et ah 
J2014h . 
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Table 2. MagAO Astrometry on 2013 April 6 (UT) 


Filter Plate Scale ("/pixel) 
r' 0.007917 ±0.000015 

i' 0.007907 ±0.000015 

± 0.007911 ±0.000012 

Ys 0.007906 ± 0.000014 


ame 

Separation (") 

P.A. 


’) 

B 

2.717 

± 

0.030 

298.8 

± 

0.6 

cc2 

1.962 

± 

0.004 

69.8 

± 

0.3 

B 

2.671 

± 

0.006 

300.0 

± 

0.3 

cc2 

1.962 

± 

0.004 

69.7 

± 

0.3 

B 

2.679 

± 

0.004 

300.0 

± 

0.3 

cc2 

1.995 

± 

0.003 

69.7 

± 

0.3 

B 

2.684 

± 

0.005 

299.9 

± 

0.3 

cc2 

1.991 

± 

0.005 

69.6 

± 

0.3 
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Fig. 1.— CT Cha in MagAO filters. Image contrasts are adjusted to bring out the objects. 
Top row: unsaturated dataset showing the PSF. Middle row: reduced, saturated dataset 
before any halo subtraction. Bottom row: primary’s halo removed by subtracting a rota- 
tionally symmetric PSF. The r' image was further smoothed to enhance the signal-to-noise 
ratio. We £t the bottom row images with the DAOPHOT allstar PSF fitting photometry 
task. At i', P, and Ys PSF fits using unsaturated images of A as the PSF of B were excellent, 
with very small (< 10%) PSF fitting errors. At r' where B is very faint (r' ~ 22 mag) the 
S/N was lower, so the fitting error increased to > 20%. 
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Fig. 2.— Seven-year relative astrometry for CT Cha B (top) and cc2 (bo ttom). Proper 


motio n of CT Cha A (/i„cos5 = —21.3 mas/yr, /i^ = 6.3 mas/yr) is from ISchmidt et al. 
(120081 ). The ellipse in dashed line shows the convolution between the orbital motion and 
astrometric uncertainty. Apparently cc2 is not a bounded nor a non-moving background 
object. It is, in fact, a background star moving northwestward at ~ 15.4 mas/yr (/!„ cos 5 ~ 
—8.2 mas/yr and ~ 13.0 mas/yr; absolute proper motion). For CT Cha B, the MagAO 
value is consistent with those measured from archive NACO images, showing that it is a 
physical companion. 
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Fig. 3.— Variation of the H 2 O-K 2 index with snrface gravity for 2000 — 2800 K BT-Settl 
synthetic spectra. We found that for log g = 3.0 to 4.5 the models give similar indices. 
Hence this is a reasonably gravity insensitive T^jf index. CT Cha B ’s index is sirnilar t o 
that of other young late M-dwarfs in the Taurus star-forming region flMuench et al.l 120071) . 
and is more consistent with the 2500 K models. We conservatively pick 2500 ± 100 K for 


CT Cha B. 
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Fig. 4.— Reduced as a function of extinction for ±0.3 mag normalization uncertainty 
and different surface gravity. Note that the degree of freedom (htting photometry at F, z', 
Ys, J, and H) is 3 because we normalized the model at Ks- If we £t the minima of these 
curves, then this analysis favors Ay = 3.4 ±1.1 mag, independent of surface gravity. The 
relatively high xl is likely due to some systematic errors that are discussed in section 3.3.2 
of the text. 
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Fig. 5.— SED fitt ing with the reddened 2500 K BT-Settl spectrum. The JHKg photometr y 
was adopted from Schmidt et ah fj2008 '). and the extinction law was from Fitzpatrick ( 1999 ). 
We matched the reddened model with the observed Ks flux, and calculated the reduced 
to determine the extinction (see Figure 0]). We only utilized the apparent fluxes in order 
to avoid the 10 — 20% uncertainty in CT Cha’s distance. The r' photometry is very much 
higher than photospheric due to Ha accretion luminosity. Note that r' was not used in fit 
of the Xr to emission. 
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Fig. 6.— Color-Magnitude diagram with DUSTY evolutionary tracks. The shaded area is 
the only area consistent with the observed Tg// (2400-2600 K), luminosity {\.og{Li^oilL q) = 
—2.68 ± 0.25), and estimated age of CT Cha B (1- 5 Myr). This yields a mass estimate of 
~ 14 — 24 Mj, consistent with lSchmidt et al.l fl2008[) . CT Cha B is very unlikely to be under 
the deuterium-burning limit, so it is a very low-mass brown dwarf. 








